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1.0  INTRODUCTION 


This  study  was  performed  to  determine  natural  frequencies 
and  mode  shapes  of  a  small  test  horn  using  a  computer  program 
COMIC  (COMplex  Acoustic  Pressure  Analysis),  which  was  developed 
by  battelle  Columbus  Laboratories  of  Columbus,  Ohio. 

COMIC  is  an  acoustic  finite  element  code  which  performs  a 
three-dimensional  steady-state  pressure  analysis  due  to  harmonic 
boundary  conditions  associated  with  an  acoustic  medium,  and 
determines  acoustic  pressure,  natural  frequencies,  and  the 
associated  mode  shapes  where  the  pressure  and  its  derivatives 
are  normalized  to  unity. 

Reference  1  contains  features  and  capabilities  of  COMIC,  and 
a  brief  statement  of  theoretical  development  and  several  example 
problems  with  listings  and  computer  printouts. 

The  full  capabilities  and  accuracy  of  the  program  were 
evaluated  in  Reference  2  by  comparing  the  finite  element  results 
to  the  analytical,  experimental,  and  other  numerical  techniques. 
Reference  2  shows  that  COMIC  provides  excellent  results  for 
problems  determining  frequencies  and  mode  shapes. 
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2.0  MODEL  FOR  ANALYSIS 

Figure  1  is  the  side  view  of  a  small-scale  facility  test 
horn  which  has  a  rectangular  cross  section  up  to  L=l80  inches 
measured  from  the  back  wall  and  has  a  square  cross  section  from 
L=  180  Inches  to  L=448  inches.  Starting  at  L=448  inches,  the  test 
horn  transitions  from  a  square  cross  section  to  a  circle  for  one 
air  modulator  configuration.  For  two  air  modulators,  the  test 
horn  divides  into  two  identical  tapered  parts,  each  45  inches 
long.  The  distance  between  the  geometric  centers  of  each 
separated  tapered  part  is  18  inches. 


2.1  First  Model  —  L-384  Inches 

Shown  in  Figure  2  is  the  finite  element  model  for  analysis. 
Two  types  of  models  were  run  for  comparison.  The  coarse  model 
has  375  grid  points  and  224  solid  elements,  while  the  finer  one 
has  620  grid  points  and  400  solid  elements.  The  results  are 
shown  in  the  first  two  columns  of  Table  1.  A  comparison  of  the 
results  indicates  that  the  difference  in  frequencies  are  negligi¬ 
ble.  Therefore,  the  coarse  model  was  used  in  running  the  three 
following  cases: 

(1)  Straight  Blocks  and  Missile:  A  straight  block  wall  17 
inches  thick  is  placed  at  20  inches  parallel  from  the  back  wall, 
reaching  half  the  height  of  the  ceiling.  In  addition,  a  95-inch- 
long  missile  is  also  placed  axially  at  X=l68  inches. 

(2)  Irregular  Blocks  and  Missile:  Blocks  are  laid  such 
that  the  back  of  the  block  wall  is  parallel  to  the  back  wall, 
while  the  front  side  is  jagged. 

(3)  Tilted  Blocks:  Blocks  are  piled  to  make  an  angle  with 
the  back  wall.  This  case  does  not  have  a  missile. 
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2.2  Second  Model  —  L=448  Inches 

This  model  is  constructed  by  adding  a  cavity  64  inches  long 
with  a  constant  cross  section  12  inches  by  12  inches  to  the 
right-end  of  the  first  model.  Two  cases  are  considered  -  no 
cavity  and  tilted  blocks  only. 

2.3  Pull  Model  —  L=503  Inches 

Figure  3  is  the  finite  element  model  for  analysis  where  a  Y 
shaped  and  tapered  section  is  included.  The  source  of  sound  is 
at  the  end  of  the  model  (X-503  inches). 


3.0  RESULTS 

3. 1  Natural  Frequencies 

Table  1  contains  the  acoustic  natural  frequencies  of  all 
models  considered  in  this  report.  The  results  from  the  coarse 
model  show  excellent  agreement  with  those  from  the  finer  model. 

We  found  that  placing  blocks  and  a  missile  inside  the  model  did 
not  make  a  significant  change  in  frequencies.  This  is  because 
the  longitudinal  mode  becomes  dominant  at  lower  frequencies  and 
the  wave  length  becomes  much  larger  than  the  thickness  of  the 
block  wall.  Another  reason  may  include  the  fact  that  the  blocks 
are  piled  up  only  half  the  height  to  the  ceiling,  and  thus  do  not 
contribute  much  in  changing  the  size  or  shape  of  the  model. 
However,  the  increase  of  the  length  of  the  model  reduces  natural 
frequencies  as  expected. 


3.2  Mode  Shape 

For  the  natural  frequency  analysis,  the  output  of  a  COMIC 
run  contains  acoustic  pressure  at  grid  points  and  its  partial 
derivatives,  and  they  are  normalized  to  unity. 

As  there  was  a  need  for  graphical  display  of  the  structural 
model  and  computer-generated  output,  a  program  was  written  to 
convert  the  COMIC  output  data  into  the  one  which  operates  on  CADS 
(Computer  Aided  Design  System)  using  DI-3000  graphics  package. 
Heference  3  presents  detailed  user  instructions  for  CADS. 

The  numerical  values  of  acoustic  pressure  differences  are 
represented  by  eigenvectors  in  CADS,  and  CADS  displays  the  eigen¬ 
vectors  mode  shape  data  at  the  node. 

From  Figures  4  through  19,  the  top  view  of  the  original 
geometry  was  drawn  using  a  solid  line,  and  mode  shapes  are  repre¬ 
sented  by  a  dotted  line  which  connects  tips  of  eigenvectors  at 
grid  points.  A  positive  pressure  difference  may  be  identified  by 
an  arrow  going  in  the  right  and  upward  direction  when  two  corre¬ 
sponding  grids  are  connected. 


A  node  point  is  the  crossing  point  of  a  solid  line  with  a 
corresponding  dotted  line,  and  a  node  line  is  constructed  by 
simply  connecting  these  node  points.  This  is  shown  in  Figures  5 
and  6.  Mode  shapes  of  the  three  models  are  plotted  in  Figures  4 
through  19  using  a  top  plate  with  an  open  window  as  an  example. 

Figure  4  shows  an  original  finite  element  model  of  the  top 
plate  in  solid  lines  and  a  mode  shape  of  the  first  model 
(L=384  inches)  at  the  first  frequency  (13.802  Hz)  in  dotted 
lines.  The  window  is  represented  by  the  area  crossed,  and  the 
pressure  difference  around  the  window  becomes  zero  as  shown  in 
the  figure. 

Determination  of  the  correct  sign  of  an  eigenvector  is 
illustrated  by  taking  as  an  example  a  point  at  the  lower  left 
corner  in  the  figure.  This  point  on  the  solid  line  is  shifted  to 
the  upper  right  position  on  the  dotted  line,  and  thus  a  vector 
connecting  these  two  points  is  positive,  and  a  positive  acoustic 
pressure  is  present  at  that  point.  As  this  is  true  over  the 
surface  of  the  plate  and  no  line  crosses  its  original  solid  line, 
the  pressure  becomes  positive  all  over  in  the  cavity,  thus  pro¬ 
ducing  mode  (0,0,0),  which  is  believed  to  be  analogous  to  a  rigid 

body  motion  in  a  structure  analysis. 

Figure  5  shows  a  mode  shape  at  the  second  frequency 
(30.637  Hz)  and  is  different  from  the  first  mode.  Following  the 
bottom  line  shows  that  the  acoustic  pressure  difference  starts 
with  a  negative  value,  but  becomes  positive  after  the  dotted  line 
crossed  over  the  solid  line,  producing  a  node  point  at 
X  =  1 5 1  Inches.  Other  node  points  are  all  at  the  same  X-coordinate 
but  at  different  Y-coordlnates .  Therefore,  a  node  line  parallel 
to  the  Y-axis  is  drawn  as  shown  in  Figure  5.  As  this  is  true  in 
all  other  X-Y  planes  along  Z-axls,  a  node  plane  can  be  identified 

as  a  Y-Z  plane  at  X=lbl.04  Inches  and  the  node  becomes  (1,0,0). 

The  third  mode  In  Figure  6  has  two  node  planes,  X=106.44 
inches,  and  X=275.8  inches.  A  mixed  node  (2,1,0)  is  shown  in 
Figure  H. 
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Mode  shapes  in  Figures  9  through  13  are  drawn  for  the  second 
model  (L-448  inches)  using  the  same  top  plate  with  an  open  window 
but  extended  longitudinally,  and  those  ±n  Figures  14  through  19 
are  for  the  full  model  (L-503  inches). 

A  comparison  of  Figures  5,  10,  and  15  reveals  that  the  first 
node  plane  moves  toward  the  direction  of  adding  elements.  In 
fact,  the  first  node  plane  is  present  at  X-151.04  Inches  for  the 
first  model  (L=384  inches)  in  Figure  5,  and  at  X-2 44.50  inches 
for  the  second  model  (L-448  inches)  in  Figure  10,  and  at  X-244.50 
inches  for  the  full  model  (L-503  inches)  in  Figure  15. 


4.0  CONCLUSIONS 

A  small-scale  facility  test  horn  was  analyzed  by  making  use 
of  the  three-dimensional  acoustic  finite  element  computer  code 
COMIC.  Acoustic  natural  frequencies  and  mode  shapes  are  obtained 
for  three  different  models. 

Mode  shapes  are  presented  graphically  using  CADS,  and  node 
planes  are  identified. 

We  found  that  placing  a  wall  of  blocks  and  a  missile  inside 
the  model  does  not  make  any  significant  change  in  natural  fre¬ 
quencies,  and  that  the  node  planes  have  moved  in  the  direction  of 
added  elements. 
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